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\-J ' 1 Introduction degree and pk is the fraction of nodes with degree k, the 

(N ■ 

,_^ ■ network having 7 < 3 is highly robust against the random 

"> One of the most important properties of complex net- . .. . , .. , ... 
, £_< ■ failure: you have to remove almost all nodes to disinte- 
rs ■ works, such as the WWW, the Internet, social and biolog- , , TT r,-^ -, 

%4 , grate the network. However, bb networks are fragile to 

C^ 1 

ical networks, is robustness against failures or intentional . . . 

the intentional attack: the network is destroyed if a small 

fraction of hubs are removed. Later, this observation has 
been analytically supported [51IWTU]. The robustness of 
networks against other percolation-like processes such as 
the betweenness-based attack [TT] and degree- weighted at- 
tacks [T5] have also been studied. 
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Abstract. We investigate robustness of correlated networks against propagating attacks modeled by a 
susceptible-infected-removed model. By Monte-Carlo simulations, we numerically determine the first crit- 
ical infection rate, above which a global outbreak of disease occurs, and the second critical infection rate, 
above which disease disintegrates the network. Our result shows that correlated networks are robust com- 
pared to the uncorrelated ones, regardless of whether they are assortative or disassortative, when a fraction 
of infected nodes in an initial state is not too large. For large initial fraction, disassortative network be- 
comes fragile while assortative network holds robustness. This behavior is related to the layered network 
structure inevitably generated by a rewiring procedure we adopt to realize correlated networks. 

PACS. XX.XX.XX No PACS code given 



attacks [1,2,3,4,5.6]. Albert et al. numerically studied 
the robustness of networks against random failure, where 
nodes are sequentially removed with equal probability, and 
intentional attack, where hubs (nodes with large degrees) 
are preferentially removed [7|. When a network has the 
scale- free (SF) degree distribution pk oc fc~ 7 , where k is 
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Other attacks to networks may occur as propagating be disassortative [T91I201 I2] . Degree correlations affect dy- 

processes such as computer viruses do [T5irm[T5] . Let namics on networks. Previous studies have reported the 

us consider the susccptible-infectcd-rcmoved (SIR) model effects of the degree correlations on various dynamics, 

with infection rate A (and recovery rate \i = 1) on a net- e.g., percolation [2T1I22] . target attack [23ll24j . Ising model 

work. The system has, in general, two critical infection [25j . susceptible- infected-susceptible model |26| . informa- 

rates A c i and A C 2 [TB1I17] . Above the first critical infcc- tion spreading [57], etc. As for the SIR model, Vazquez 

tion rate A c i, a global outbreak of disease occurs. In a and Moreno showed that the increase of the assortativity 

global outbreak, a fraction of the nodes become infected makes easier to induce a global outbreak J2E|. However, 

and eventually removed. The remaining network of susccp- we do not know the effect of the degree correlation of nct- 

tible nodes, however, may survive as a giant component works on the robustness against such propagating attack 

up to the second critical infection rate A C 2, above which modeled by the SIR model. 

the network is finally disintegrated. Thus we here adopt In this paper, we perform a Monte-Carlo simulation for 

A C 2 as a measure of robustness against propagating at- the SIR model on networks having a degree correlation. 

tacks |17j . Note that A C 2 is in principle equal to or larger We numerically determine A c i and A C 2 to show that cor- 

than A c i. related network is robust compared to the uncorrelated 

t rroi i-i. c i. i • c i.- \ • one regardless of whether it is assortative or disassorta- 

In [18J, the first critical infection rate A c i is approxi- ° 

, i j . i t i • r i- tive, when a fraction of infected nodes in an initial state 
mately derived. In particular, any positive infection rate 

i , , , ,, , . n , , , is not too large. For large initial fraction, disassortative 

induces a global outbreak in Si 1 networks having 7 < 3, b ° ' 

., n at j- j ii nm j 1 network becomes fragile while assortative network holds 

meaning A c i = U. iNcwman studied the MR model on un- ° 

, , , , . P , . .■,.■,., 1 robustness. This behavior is related to the layered network 

correlated networks m terms of transmissibinty to snow J 

, , , , ^ n c ^o nvn tt 1 »/r j structure inevitably generated by a rewiring procedure we 

that A C 2 > even for 7 < 3 [16J . Hasegawa and Masuda ° J ° 

, 1 , , . , . , , , ,, 1, x . adopt to realize correlated networks l25ll29l . 

extended his analysis to evaluate the robustness A C 2 m ^ ' — " — ' 



some vaccinated networks [17] . 

Real networks often have some degree correlations [THl 



2 Model 



[20] . A network is said to be assortative when nodes with We consider the SIR model on a given network with TV 

similar degrees tend to connect with each other. A network nodes. Each node takes one of the three states: suscep- 

is called disassortative when nodes with high degrees tend tible, infected, or removed. A fraction p soo d of the nodes 

to connect to nodes with low degrees. Empirical data in- are initially infected and other nodes are susceptible. Sus- 

dicates that social networks are likely to be assortative, ceptible nodes get infected at a rate proportional to the 

while technological networks and biological networks to number of infected neighbors: the susceptible node gets in- 
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fected with probability XAt within a short time At when An uncorrclatcd network having the above distribu- 
it is adjacent to an infected node. An infected node be- tions can be generated by the configuration model [2J. 
comes removed, irrespective of the neighbors' states, at a The network so obtained is then randomly rewired so as 
unit rate \x = 1, i.e., with probability At within short time to have a correlation without changing the degree distri- 
ct, bution. We adopt the following scheme proposed in [25] 
In the final state, each node takes either susceptible or [25] : (a) Select randomly two edges and look up the de- 
removed state. We call the connected components of re- grees of four nodes linked by these edges, (b) Rewire the 
moved nodes and susceptible nodes the infected networks two edges, with probability p r0 wirc, in a way that one edge 
and the residual networks, respectively. Then, the critical links two nodes having the higher degrees and the other 
infection rates A c i and A C 2 are given in terms of percola- links the remaining two nodes, and otherwise rewire them 
tion transition. We define R max (N) and S max (N) to be the randomly, (c) Repeat (a) and (b) until the degree corre- 
mean largest component sizes of the infected networks and lation (see below) reaches a stationary value. To generate 
residual networks, respectively. Noting that an infection a disassortative network, we replace the step (b) with the 
quickly dies out (spreads globally) for A < A c i(> A c i) and following step: (b)' Rewire the two edges, with probability 
the residual networks include a giant component (consist p 10 wirc, in a way that one edge links two nodes having the 
of only finite components) for A < A c2 (> A C 2), we expect highest and lowest degrees and the other links the remain- 
that Rmax(N) and S maiK (N) behave as ing two nodes, and otherwise rewire them randomly. 

{const., iv} (A < A c i) Here we adopt the following quantity for a measure of 

{.Rmax, S max } oc ^ {N,N} (A c i < A < A c2 ) ■ (!) the degree correlation [25]: 

{TV, const.} (A > A c2 ) 

In the following sections, we perform Monte-Carlo sim- , _ z^fc £** — Z-ifc ^kk /a\ 

1 2-ik °kk 

ulations for the SIR model on the uncorrclatcd and cor- 
related networks. Here we adopt the following degree dis- 

where tkk is the probability that both ends ol a randomly 

tributions: an exponential degree distribution, 

chosen edge have degree fc, and t kk is that of the uncorrc- 

1 /2\ fc ~ 2 
Pk = - [ - ) , (2) lated networks (with the same degree distribution) . Uncor- 

rected networks take A = 0, and a network is regarded 
and a SF degree distribution with 7 = 3, 

as assortativc (disassortative) when A > 0(< 0). After 
Pk = 12 oc fc~3_ ^ 

k(k + l)(fc + 2) the above process, the mean value of A takes a p r cwirc- 

Note that those networks have the same average degree dependent value. Choosing p rC wirc appropriately, we have 

(fc) =4. a network with a desirable degree correlation. 
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Fig. 1. Results for uncorrelated networks with an exponen- 



3 Result: Uncorrelated Case 

In this section, we show the results of Monte-Carlo simu- 
lations for the SIR model on uncorrelated networks. The 
number of graph realizations is 500 and the number of SIR 
runs on each realization is 100. The number of nodes is 
taken from N = 2 9 to 2 14 . 

First, we consider uncorrelated networks with the ex- 
ponential degree distribution @. The order parameters 
Rm&x(N)/N and S max (N)/N with several network sizes 
are shown in Fig. [TJa). Here we set p scc d = 0.1. As the 
infection rate A increases, R ma , x {N) / N monotonically in- 
creases while S max (N)/N monotonically decreases. 

To determine the critical points, we introduce the frac- 
tal exponents ipn and ips of the infected and residual net- 
works defined by [5U1I5T] 

SmaxW cx A^ S , R max (N) cc N** , (5) 

respectively. Noting that Eq.([T]), we expect 

{0,1} (A<A c i) 

{Vfl, c , 1} (A = A cl ) 
{ip R ,ips}<x\ {1,1} (A cl <A<A c2 )- (6) 

{l,i>s,c} (A = A c2 ) 
{1,0} (A>A c2 ) 

For numerical computations, we evaluate the effective ex- 
ponents of finite N: 



dlogiV ™ v ' dlogJV 

tial degree distribution: (a) order parameters R m ax(N)/N and mi i n\r\ j / i i\r\ -ii i • + \ 

& v ; H v ;/ Then, ipR[N) and ips{N) with several sizes cross at A 

Sma,x(N)/N, (b) fractal exponents ips(N) and tpn(N), (c) scal- 
ing plot of Rmax(N) and (d) scaling plot of S l ma x(A r ). The num- 
ber of nodes iV is 2048 (red-circle), 4096 (green-triangle), 8192 
(blue-diamond), and 16384 (purple-square). The full and open 
symbols denote the results for the residual and infected net- 
works, respectively. 



A c i and A = A c2 , respectively. 

By using data of Fig. [lja) , we plot the effective fractal 
exponent in Fig.QJb). We have a crossing point of ipn(N) 
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at A c i « 0.05, and of i/js(N) at A C 2 ~ 0.72. At the crossing 
points, the corresponding fractal exponents take the same 
value ipR.c — ips.c = 2/3. The property that the largest 
component size at criticality is proportional to N 2 ' 3 is 
also observed in the mean field SIR model [521155] and 
percolation of the random graph [34] . 

We further determine the critical exponent by using a 
finite size scaling. For R max (N), we assume the following 
scaling form [31] : 



(a) 



where f{x) is a scaling function satisfying 
r x i-4>r,c f or 3; > l 



(8) 



(9) 



I const for x <C 1 

AX = A — A c i, and (3 is the critical exponent of the or- 
der parameter, Rm ax (N)/N oc |^Ap. The scaling form 
for S' max (A r ) is assumed in a similar way. Our scaling re- 
sults for Rmax(N) and S max (N) are shown in Figs. QJc) 
and (d). For both plots, the fitting parameter /3 is taken 
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as (3 = 1. This indicates that these transitions are in the Fig. 2. Results for uncorrelated SF networks: (a) order pa- 
mean field universality class. Similar scalings for the mean rameters Rma.*{N)/N and SmiLx(N)/N, (b) fractal exponents 
cluster sizes of the residual and infected networks are also 4>s(N) and i/jr(N), and (c) scaling plot of 5 max (A/'). The num- 
assumed in a similar way [3T]- We have a good data col- ber of nodes N is 2048 (red-circle), 4096 (green-triangle), 8192 
lapse for both by setting /3 = 1 (not shown). (blue-diamond), and 16384 (purple-square). The full and open 

The results for the SF network are shown in FigH As s y mbols denote the results for the residual and infected net " 
N increases, i()r(N) near A = approaches one, indicating ' 

that A c i = for SF networks with 7 = 3 [TF . On the other 
hand, ips(N) with several sizes cross at A C 2 ~ 0.56, which 
is smaller than that of the exponential distribution. This 
is consistent with the fact that the SF network with 7 = 3 
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is fragile compared to the random graph with the same 
average degree |17j . The same analysis is performed with 
several initial seed fractions p scc d — 0.3, 0.5, 0.7. Our simu- 
lations indicate that the SF network remains fragile in this 
range of p SO od- In FigHJc), our scaling result shows that the 
transition of residual networks belongs to the mean field 
universality class, /3 = 1, even when the network is scale- 
free with 7 = 3. Here we assumed 4>s,c = 2/3. Although 
our numerical data (Fig|2Jb)) seems to be less precise due 
to small samples, we performed a similar analysis for a 
dynamic percolation in the same universality as the SIR 
model to show precisely ips,c = 2/3 for that model (not 
shown) . 

4 Result: Correlated Case 

We performed similar simulations on correlated networks 
with several degree correlation A. Here A takes a value in 
a range [—0.1,0.5]. 

The critical infection rates A c i and A C 2 of the net- 
work with the exponential degree distribution are plotted 
as a function of A (FigJ3]). As expected in [2S], A c i de- 
creases when A increases (Figj3]Ja)). On the other hand, 
we observe a nontrivial behavior for A C 2 (Fig(3Jb)). For 
Psoed = 0.1, 0.3, 0.5, the uncorrelated network is most frag- 
ile: A C 2 increases when the network is correlated, regard- 
less of whether it is assortative or disassortative. This ten- 
dency decreases with the increase of p se cd, and we did not 
observe such a behavior for p soc d = 0.7. There, a disassor- 
tative network is fragile than the uncorrelated one. Figure 
0] is the result for SF network. We can observe a similar 
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Fig. 3. (a) Aci and (b) A C 2 of the correlated network with 
an exponential degree distribution as a function of degree 
correlation A. The initial seed fraction p S ocd is set to 0.1 
(red-circle), 0.3 (green-triangle), 0.5 (blue-diamond), and 0.7 
(purple- square) . 

dependence of A C 2 on A. Our scaling results indicate that 
^S,c (and (3) tends to decrease (increase) with the increase 
of the assortativity, although it was hard for our simula- 
tion to determine precise values because of large error-bars 
(not shown). More extensive simulation or another effec- 
tive method should be performed for the evaluation of the 
critical properties in detail. 

The above behavior of A C 2 is due to the layered network 
structure of correlated networks [25] . In Fig. [SJ we profile 
the adjacency matrix of SF network for several A. Assorta- 
tive networks consist of monolayers of nodes where nodes 
with similar degrees are connected to each other, while dis- 
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n 5 Summary 




Fig. 4. A C 2 of the correlated network with the SF degree 
distribution as a function of degree correlation A. The initial 
seed fraction p see d is set to 0.1 (red-circle), 0.3 (green-triangle), 
0.5 (blue-diamond), and 0.7 (purple-square). 

assortative networks consist of bilaycrs where nodes with 
the low-degree side are connected to a range of high-degree 
nodes. For such a layered structure, there are few edges 
between different layers. 

Our propagating attack consists of two elements: the 
random initial failure realized by p se ed and the SIR dy- 
namics starting from infected nodes. When the SIR dy- 
namics starts from a single infected node, the outbreak 
from the node, local or global, is likely to be confined in 
the layer containing the node and all other layers will re- 
main as the giant component after the attack. Therefore, 
correlated networks, assortative or disassortative, are ro- 
bust compared to uncorrelated ones for small p see d- For 
large p SC od, on the other hand, the initial failure damages 
a number of susceptive nodes and the remaining network 
becomes marginally percolating at the beginning. In this 
situation the more assortative (disassortative) is the cor- 
relation, the more robust (fragile) is the network against 
the random failure, just as seen in the percolation on cor- 
related networks [25] . 



We have investigated the robustness of correlated net- 
works against a propagating attack modeled by the SIR 
model. Our numerical results show that correlated net- 
works are robust compared to the uncorrelated ones, re- 
gardless of whether they are assortative or disassortative 
when a fraction of infected nodes in an initial state is not 
too large. This behavior is related to the layered struc- 
ture of our correlated networks. For large p SC cd, disassor- 
tative network becomes fragile. Note that the layered net- 
work structure is not an intrinsic property of correlated 
networks (e.g., we can construct a disassortative network 
having almost no layered structure |22j). It remains open 
for future works to search a mechanism of strengthening 
the robustness of real correlated networks against propa- 
gating attacks. 

Previous studies for epidemics on complex networks 
have mainly focused on the (first) critical infection rate 
A c i, above which a global outbreak occurs. Instead, we 
have focused on the second critical infection rate A C 2 in 
this paper. In particular, it is important when we con- 
sider epidemics on community networks, where we should 
distinguish local epidemics, i.e., an outbreak confined in 
a single community, and global epidemics, i.e., epidemic 
spreading through the whole communities |35j . We expect 
that the robustness against propagating attacks will be a 
good measure for searching effective vaccines in commu- 
nity networks. 
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Fig. 5. The ordered adjacency matrix A of a SF network with (a) A — —0.1, (b) A — 0, and (c) .4 ~ 0.8. The number of 
nodes is N = 1000. The matrices have N x N entries Aij, where Aij = 1 (black) if nodes i and j are connected and Aij — 
(white) otherwise. The nodes are arranged in the order of increasing degree. 
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